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Aim 



The aim of the percolation theory course is to provide a challenging and stimulating introduction 
to a selection of topics within modern theoretical condensed matter physics. 

Percolation theory is the simplest model displaying a phase transition. The analytic solutions 
to Id and mean-field percolation are presented. While percolation cannot be solved exactly for 
intermediate dimensions, the model enables the reader to become familiar with important concepts 
such as fractals, scaling, and renormalisation group theory in a very intuitive way. 

The text is accompanied by exercises with solutions and visual interactive simulations for the 
percolation theory model to allow the readers to experience the behaviour, in the spirit "seeing is be- 
lieving". The animations can be downloaded via the URL http:/ /www. cmth.ph.ic.ac.uk/kim/cmth/ 

I greatly appriciate the suggestions and comments provided by Nicholas Moloney and Ole Peters 
without whom, the text would have been incomprehensible and flooded with mistakes. However, 
if you still are able to find any misprints, misspellings and mistakes in the notes, I would be very 
grateful if you would report those to k.christensen@ic.ac.uk. 
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1.1 Introduction 



Percolation theory is the simplest not exactly solved model displaying a phase transition. Often, 
the insight into the percolation theory problem facilitates the understanding of many other physical 
systems. Moreover, the concept of fractals, which is intimately related to the percolation theory 
problem, is of general interest as it pops up more or less everywhere in Nature. The knowledge of 
percolation, fractals, and scaling are of immense importance theoretically in such diverse fields as 
biology, physics, and geophysics and also of practical importance in e.g. oil recovery. We will begin 
gently by developing a basic understanding of percolation theory, providing a natural introduction 
to the concept of scaling and renormalisation group theory. 

1.2 Preliminaries 

Let P(A) denote the probability for an event A and P(A% D A2) the joint probability for event A\ 
and A2. 

Definition 1 Two events A\ and A<i are independent 4^ P{A\ C\ A2) = P{Ai)P{A2). 

Definition 2 More generally, we define n > 3 events A\,A2, ■ ■ ■ ,A n to be mutually independent 
if P{A\ n A2 n • • • fl A n ) = P{Ai)P{A2) ■ ■ ■ P{A n ) and if any subcollection containing at least two 
but fewer than n events are mutually independent. 

Let each site in a lattice be occupied at random with probability p, that is, each site is occupied 
(with probability p) or empty (with probability 1 —p) independent of the status (empty or occupied) 
of any of the other sites in the lattice. We call p the occupation probability or the concentration. 

Definition 3 A cluster is a group of nearest neighbouring occupied sites. 

Percolation theory deals with the numbers and properties of the clusters formed when sites are 
occupied with probability p, see Fig. (1.1). 




Figure 1.1: Percolation in 2d square lattice of linear size L = 5. Sites are occupied with probability 
p. In the lattice above, we have one cluster of size 7, a cluster of size 3 and two clusters of size 1 
(isolated sites). 

Definition 4 The cluster number n s {p) denotes the number of s- clusters per lattice site. 

The (average) number of clusters of size s in a hypercubic lattice of linear size L is L d n s (p), d being 
the dimensionality of the lattice. Defining the cluster number per lattice site as opposed to the 
total number of s-clusters in the lattice ensures that the quantity will be independent of the lattice 
size L. 
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For finite lattices L < oo, it is intuitively clear, that if the occupation probability p is small, 
there is only a very tiny chance of having a cluster percolating between two opposite boundaries 
(i.e., in 2d, from top-to-bottom or from left-to-right). For p close to 1, we almost certainly will 
have a cluster percolating through the system. In Fig. 1.2, sites in 2d square lattices are occupied 
at random with increasing occupation probability p. The occupied sites are shown in gray while 
the sites belonging to the largest cluster are shown in black. Unoccupied sites are white. Note that 
for p ~ 0.59, a percolating cluster appears for the first time. 




Figure 1.2: Percolation in 2d square lattices with system size L x L = 150 x 150. Occupation prob- 
ability p = 0.45,0.55,0.59,0.65, and 0.75, respectively. Notice, that the largest cluster percolates 
through the lattice from top to bottom in this example when p > 0.59. 

Definition 5 The percolation threshold p c is the concentration (occupation probability) p at which 
an infinite cluster appears for the first time in an infinite lattice. 

Note, that p c is defined with respect to an infinite lattice, that is, in the limit of L — > oo. Table 
(1.1) lists the percolation threshold in various lattices and dimensions. 

Exercise 1 Why is p c not well defined in a finite lattice ? 

1.3 Percolation in Id 

We will consider the percolation problem in Id where it can be solved analytically. Many of the 
characteristic features encountered in higher dimensions are present in Id as well, if we know 
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Table 1.1: The percolation threshold for the site percolation problem is given in column 3 for various 
lattices in various dimensions. Column 2 lists the number of nearest-neighbours (nn), also known 
as the coordination number. Within a given dimension, the percolation threshold decrease with 
increasing number of nearest-neighbours. The site percolation problem has a counterpart called 
the bond percolation problem: In a lattice, each bond between neighbouring lattice sites can be 
occupied (open) with probability p and empty (closed) with probability (1 — p). A cluster is a group 
of connected occupied (open) bonds. NB: In all cluster is defined as a group of nearest 

neighbouring occupied sites (bonds). Note that the percolation threshold for the site-percolation on 
high-dimensional hypercubic lattices, where loops become irrelevant, approaches that of the Bethe 
lattice l/(z — 1), if we substitute the coordination number z with 2d. 

where and how to look. Thus the Id case serves as a transparent window into the world of phase 
transitions, scaling, scaling relations, and renormalisation group theory. 

Imagine a Id lattice with an infinite number of sites of equal spacing arranged in a line. Each 
site has a probability p of being occupied, and consequently 1 — p of being empty (not occupied). 
These are the only two states possible, see Fig. 1.3. 

x • • y, • x • • • • • x x • x « x 

Figure 1.3: Percolation in a Id lattice. Sites are occupied with probability p. The crosses are empty 
sites, the solid circles are occupied sites. In the part of the infinite Id lattice shown above, there is 
one cluster of size 5, one cluster of size 2, and three clusters of size 1. 

What is of interest to us now and in future discussions in higher dimensions is the occupation 
probability at which an infinite cluster is obtained for the first time. A percolating cluster in Id 
spans from — oo to +oo. Clearly, in Id this can only be achieved if all sites are occupied, that is, 
the percolation threshold p c = 1, as a single empty site would prevent a cluster to percolate. 

A precise "mathematical" derivation of p c = 1 in Id goes as follows. 

Definition 6 Let n(p, L) denote the probability that a lattice of linear size L percolates at concen- 
tration p. 
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Combining the two definitions (5) and (6), we have, 



lim U(p, L) 

L— >oo 



for p < p c 

1 for p > p c . 



Consider a finite Id lattice of size L where each site is occupied with probability p. As the 
events of occupying sites are independent, all sites are occupied with probability U(p, L) = p L , see 
Fig. (1.4), and 

J for p < 1 



implying p c = 1. 



lim L) = lim p — , 

L^oo L^oo 1 for p = 1, 




Occupation probability p 



Figure 1.4: The probability of a lei lattice of linear size L to percolation at occupation probability 
p. In the limit L — * oo, U(p, L) converges to a discontinous step function. 

Let us consider the clusters formed in a Id lattice. A cluster of size s is formed when s sites 
are occupied next to one another bounded by two empty sites, see Fig. 1.3. When L — * oo, we 
can ignore the effects of the boundary sites of the lattice and the probability of an arbitrary site 
(occupied or not) being, say, the left hand side (LHS) of an s-cluster is 

n s {p) = (1 - p)p s (l - p) = (1 - p) V- (1.1) 

This expression is obtained from the assumption that the occupancy of each site is independent 
of the state of any other site. If this was not the case then it would be much more complicated. 
Note that, since all sites have equal probability of being occupied (or empty), the probability that 
an arbitrary site is part of an s-cluster is s times the probability of it being the LHS of the cluster. 
We can re-write the cluster number Eq.(l.l) for Id percolation as 

n s (p) = (l-p)V 

= (1 -p) 2 exp(ln(p s )) 
= (1 — p) 2 exp(s ln(p)) 

= (p c -p) 2 exp(--) (1.2) 
H 

with the definition of a cutoff cluster size or characteristic cluster size 

H = T~T\ = Tl T vT ~* — " — = (Pc-p)~ 1 for p^pc, (1.3) 

ln(p) ln(p c -{Pc-p)) Pc-P 
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where we, to obtain the limit, have used p c = 1 and the Taylor expansion 

ln(l — x) = —x — -x 2 — -x 3 — • • • ps —x 
where the last approximation is valid for x 0, see Fig. 1.5. 




t i u ~i — i — n — — — n 

0.0 0.2 0.4 0.6 0.8 1.0 10" 10" 10"" 10" 10 

Occupation probability p p c - p 



Figure 1.5: Percolation in Id. (a) The cluster number distribution n s (p) = (p c — p) 2 exp(— j-) for 
various values of p approaching p c = 1. The vertical lines indicate the cutoff cluster size sg(p). (b) 
By plotting s 2 n s (p) versus s/s^ s(p c — p), all the data collapses onto a function x 2 exp(— x). (c) 
The characteristic cluster size = — lr ^ p ^ diverges when p — > p c = 1. (d) In the limit p —> p c = 1, 

S£ (x {pc-pY 1 - 



Exercise 2 Verify the Taylor expansion o/ln(l — x) around x = given above. 

Thus the cutoff cluster size diverges for p — > p c as a power law in the distance from the critical 
occupation probability p c , see Fig. 1.5. The divergence of the cutoff cluster size when p — > p c is also 
seen in higher dimensions where, however, another numerical exponent will describe the divergence. 
Thus it is natural to introduce a symbol for the exponent. 

Definition 7 The critical exponent a is defined by 

S£<x\p c -p\~« forp^p c . (1.4) 
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In Id percolation theory, a = 1 and log = constant — - log \p c — p\, see Fig. 1.5. 

Let us continue our journey into Id percolation theory. For p < p c we can state that the 
probability that an arbitrary site belongs to any (finite) cluster is simply the probability p of it 
being occupied. Since the probability that an arbitrary sites belongs to an s-cluster is given by 
sn s (p), we arrive at 

oo 

J2 sn s(.P)=P ioip<p c - (1.5) 

s=l 

Using the formula for summing a geometric series, we can satisfy those who prefer a rigorous 
mathematical proof: 



oo 



J2sn s (p) = J2s(l-p) 2 p s 

s=l s=l 



00 J/^S\ 



= (1-p) 2^P 



8=1 



dp 



A I 00 



vs=l 



= p- 

How large on average is a cluster or, equivalently, how large is a cluster on average to which 
an occupied site belongs? The probability that a site is occupied is p. The probability that an 
arbitrary site belongs to an s-cluster is sn s (p). Thus the probability w s that the cluster to which 
an occupied sites belongs contains s sites is 

sn s (p) sn s (p) 

W s = 



P J2s=isn s {p) 
Thus, the mean cluster size or average cluster size S(p) is given by 



s( P ) = E 



= E 



oo 

SW S 
s=l 

s 2 n s (p) 



{EZin s (p)s 

-I oo 
1 8=1 



p V dp 



\s=l 



where the operator 



d\ 2 ( d \ ( d\ , 2 d 2 
P-r) = \P-r \\P-r \tp ~n>- 



dp J \^ dp J Y' dp J ' *" dp 2 
Using the formula for summing a geometric series and the operator p-^ twice, we finally arrive at 

1-p p c -p 

where the last equality follows because in Id the critical occupation probability p c = 1. 
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Exercise 3 Derive the result in Eq.(1.6) for the mean cluster size in Id. 

We thus see that the mean cluster size diverges for p — ► p~ , where the minus sign signifies that we 
are approaching p c from below, which is what we intuitively expect if considering an infinite lattice. 
It is not possible to approach p c from above in Id as p c = 1. This is actually the main difference 
between Id and higher dimensions, where, p c < 1 and we can approach p c from both below and 
above. 

In order to investigate in detail how the mean cluster size diverges, when taking the limit 
p — > p~ , we note that the numerator in Eq.(1.6) approaches 2p Cl so 

S ip ) = E^±P^JP^K( Pc -p)- 1 farp-p". (1.7) 

Pc - P Pc-P 

Thus in Id, the mean cluster diverges like a power law in the quantity (p c — p) when p — > p c , see 
Fig. 1.6. The same phenomenon will be encountered in higher dimensions. 

Definition 8 T/ie critical exponent 7 is defined by 

S{p) oc|p c -p|~ 7 forp^p c . (1.8) 

In Id percolation theory, 7=1. 
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Figure 1.6: Percolation in Id. The average cluster size S(p) = (1 +p)/(l — p) 
diverges when p — > p c = 1. In the limit p — ► p c = 1, S'(p) — > 2/(p c — p). 
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Definition 9 TTte correlation function or pair connectivity g(r) is the probability that a site at 
position r from an occupied site belongs to the same finite cluster. 

Note this definition excludes the contribution from the infinite cluster. That need not worry us in 
Id, where all clusters are finite if p < p c = 1. Let r = |r|. Clearly, g(r = 0) = 1, since the site is 
occupied by definition. In Id, for a site at position r to be occupied and belong to the same (finite) 
cluster, this site and the (r — 1) intermediate sites must be occupied, leaving 

g(r) =p r , 

for all p, which can also be expressed in the form 



g(r) = exp(ln(p r )) = exp(rln(p)) = exp (-^j , 



(1.9) 
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where 



£ = -rrr = w ? vt^t t = (pc-p) 1 forp-»p c = i, (1.10) 

ln(p) ln(p c - (p c - p)) (p c - p) 

where we use the expansion ln(l — x) ~ — x for small x. The quantity £ is called the correlation 
length which diverges for p — ► p c . The same phenomenon will be encountered in higher dimensions. 
Note that in Id we have = £ which is why we haven't bothered displaying a figure for £, as it 
would be identical to Fig. 1.5. However, this identity will not be true in higher dimensions, where 
we will find oc £ D , where D is the fractal dimension, but more about this later. 

Definition 10 The critical exponent v is defined by 

icc\p c -p\~ u forp^p c . (1.11) 



In Id percolation theory, v = 1. 

By summing over all possible lattice sites r of the correlation function, the mean cluster size 
can be shown to be 

£</(r) =<?(p). (1.12) 

r 

In Id, this sum is straight forward identifying r with r = 0, ±1, ±2, . . ., see Problem 2, where you 
will also discover, that the sum rule is valid in all dimensions d. 

The general pattern for the exact solutions of the Id percolation problem is that certain quan- 
tities, such as the cutoff cluster size s^, the mean cluster size S(p), and the correlation length £ 
diverge at the percolation threshold. The divergence can be described by simple power laws of 
the distance from the critical occupation probability |Ap| = |p c — p\, e.g. £ oc (p c — p) _1 , at least 
asymptotically close to p c where Ap is small. The same phenomena will be encountered in higher 
dimensions even though we cannot obtain exact analytic solutions. 

1.4 Percolation in the Bethe Lattice 

The percolation problem can be solved analytically in d = 1 and d = oo. The infinitely dimensional 
case is synonymous with the Bethe lattice, a special lattice where each site has z neighbouring 
sites, such that each branch gives rise to z — 1 other branches, see Fig. 1.7. 




site 




neighbour y 

/ 0» subbranch 



branch 



Figure 1.7: The Bethe lattice with z = 3. Each site has three neighbours. Each branch contains 
z — 1 = 2 subbranches. 
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The Id case is effectively a Bethe lattice with z = 2. Why does the Bethe lattice correspond to 
the spatial dimension d = oo you might rightly ask! Well, in a hypercubic lattice, (a) the number 
of surface sites relative to the total number of sites approaches a constant when d — ► oo and (b) 
there are no closed loops when d — ► oo. The Bethe lattice has both these properties. 

(a) Let g denote the generation, that is, the distance from a "centre site". Note, however, that in 
an infinite Bethe lattice, all sites are equivalent, so the notion of a "centre site" is not to be 
taken literally. In the figure above, the first "ring" of three sites belong to generation g = 1, 
the second "ring" of six sites belongs to the second generation and so on. The total number 
of sites in a Bethe lattice consisting of g generations is 

Total no. sites = 1 + 3 • (1 + 2 + • • • + 

1 - 29 

= 1 + 3 



1 - 2 
= 3 -2 9 -2, 

while the number of surface sites is 3 • 2 ff_1 . Thus 

No. of surface sites 3 • 2 9 ~ 1 1 
Total no. of sites ~~ 3-29-2 2 

and the surface/ volume tends to a constant. 

Exercise 4 Show, that for a general Bethe lattice with coordination number z 

No. of surface sites z — 2 



for g — > oo. 



Total no. of sites 



for g — > oo. 



In a hypercubic lattice of linear size L, the surface is proportional to the volume only when 
d — > oo: the surface in d dimensions is proportional to L d_1 while the volume is proportional 
to L d leaving 



d-l i _ 1. 

Surface oc Volume <t = Volume <* , 



that is, the surface is proportional to the volume if d — > oo. 

(b) There are no closed loops in a Bethe lattice. Starting from the "centre site" going outwards, 
one will never return to the starting point. In a hypercubic lattice, the chance (probability) 
of having a loop approaches zero as the dimension d — > oo: As an example, let us place four 
particles in a chain in a hypercubic lattice with dimension d. When the first particle has 
been placed, there are 2d nearest neighbour sites, where the second particle can be placed. 
However, for the third and fourth particle, there are only 2d — 1 possible sites, implying a 
total no. of different chains 2d • {2d — l) 2 . Calculating the number of ways to place four 
particles in a loop, we arrive at 2d ■ {2d — 2) • 1, that is, the probability of having a four loop 

No. loops 2d ■ {2d - 2) • 1 {2d - 2) n r , 

- for d —> oo. 



Total no. chains 2d ■ {2d - l) 2 {2d - l) 2 

For d = oo, there is no chance of having a closed loop, which is intuitively clear, isn't it? 

What is the critical occupation probability in a Bethe lattice, that is, at which occupation 
probability does an infinite cluster (path) appear for the first time in an infinite Bethe lattice? 
Starting from a "centre site" and going outwards, we encounter {z — 1) new neighbours. Thus, 
on average, we have p{z — 1) new occupied sites on which we can continue the path. The critical 
occupation probability is determined by the equation 

p c {z - 1) = 1 <^ p c = — !—. 

Z — I 
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Notice, that for z = 2, the Id percolation problem, we recover p c = 1. For z = 3, e.g., p c = ^. Note 
that in a hypercubic lattice with d — > oo one would expect p c — > l/(2d — 1) in order to be able to 
continue walking along a path of occupied sites. 

Definition 11 The strength of the infinite cluster P(p) is the probability that an arbitrary site 
belongs to the infinite cluster. 

Exercise 5 Discuss why, obviously, P(p) = for p < p c in an infinite lattice. 

In order to calculate P(p) for p > p c in the Bethe lattice, we introduce the quantity Q as the 
probability that an arbitrary site is NOT connected to infinity through a fixed branch originating 
at this site. Restricting ourselves to the lattice with z = 3 and using basic probability theory, we 
have that the strength P{p) of the infinite network, that is, the probability that an arbitraryly 
selected site is connected to infinity by occupied sites is 

Pip) = (Prob. site is occupied) x (Prob. at least ONE branch lead to infinity) 

= p(i-Q 3 )- 

Next, we need to determine Q, the probability that a branch originating at that site does not 
lead to infinity. This depends on whether or not the neighbouring site is occupied. We find 

Q = (Prob. site is empty) + (Prob. site is occupied) x (Prob. no subbranch leads to infinity) 

= (i-p)+pQ 2 - 

We have relied on the fact that all sites in a Bethe lattice are equivalent, so Q is also the probability 
that a subbranch is not connected to infinity. The quadratic equation can easily be solved: 



l±V(2p-l) 2 J 1 
Q = % = \ 

Exercise 6 Show that, in a general Bethe lattice where each site has z neighbours (Problem 4), 

r i 

Q = { -I 2p(z-l)-2 
{ 1 " p(*-l)(*-2)- 

When p < p c , there are no infinite clusters, by definition, so with probability 1 there is no connection 
to infinity. Thus Q = 1 is the trivial solution associated with P(p) = for p < p c . The non-trivial 
solution Q = is associated with p > p c , leaving 

!0 for p < p c 

p(l-(V0 3 ) for p>p c . 

Using a Taylor expansion for P(p) around p = p c = \, it can be shown that (see Problem 4) 

P(j>) oc [p — p c ) for p — > p+ . 
Definition 12 The critical exponent (3 is defined by 

P{P) oc (p - Pcf forp^p+. (1.13) 
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Figure 1.8: The strength P(p) of the infinite cluster in a Bethe lattice with z = 3. When p — * p+, 
the strength P(p) — > 6(p — p c ) implying (3 = 1 in a Bethe lattice. 



In the Bethe lattice, = 1, see Fig. 1.8. 

The strength P{p) measures how large a fraction of the sites in the lattice belongs to the infinite 
cluster (in an infinite lattice) and is known as the order parameter. The phenomenon that the order 
parameter (strength) becomes non-zero for p > p c is known as a phase transition and p = p c is 
often called a critical point or critical occupation probability. 

Exercise 7 What is the strength P{p) in the Id percolation problem? 

What is the mean cluster size S(p) to which an occupied site (e.g. the "centre site") belongs? Let 
T be the mean cluster size in one branch. 

S{p) = average cluster size to which the origin belongs 

= (contribution from origin) + (contributions from the 3 branches) 
= 1 + 3T. 

Thus, we have to determine T. Again, as all sites are equivalent in an infinite Bethe lattice, T will 
also denote the mean cluster size in the subbranch, so we can argue as follows: If the neighbour 
to the "centre site" is empty, probability (1 — p), there is no contribution to T from this branch. If 
the neighbour to the "centre site" is occupied (probability p) , it contributes its own mass (unity) 
to the cluster, and adds the mass T for each of its z — 1 = 2 subbranches. Thus, 

T = (1 - p) • + p • (1 + 2T) 4» T = - P forp<p c . 



l-2p 



Therefore, the mean cluster size S(p) is given by 



= 1 + 3T = i% = = 25^5) = r2fe - pr '• <li4) 

where T2 = is the amplitude, and the critical exponent 7 = 1 in the Bethe lattice. This is the 
exact result for a mean cluster of size S(p) and we notice that it diverges for p — ► p c . 

An occupied site either belongs to the infinite cluster or a finite cluster. Thus, the generalisation 
of Eq.(1.5) is valid for all p 



00 



P(p) + sn s (p) = p Vp. (1-15) 



13 



The sum runs over all finite s and excludes the infinite cluster. Of course, Eq.(1.15) reduces to 
Eq.(1.5) for p < p c since the strength P(p) is identical zero in this range. 

Exercise 8 Reflect upon the identity Eq.(1.15) for p = p c . What is the sum J2^=i sn s {p c ) ? Is your 
answer also valid in d = 1? 

Critical phenomena occur in a variety of systems. The Bethe lattice approximation for percolation 
theory is somewhat analogous to the mean-field approximation for magnetic systems dislaying a 
phase transition, or the Van der Waals equation for a fluid. In all cases, the order parameter 
becomes non-zero above (or below) the critical point. In a magnetic system, the order parameter 
is the spontaneous magnetisation per spin m, which is zero above the critical temperature (Curie 
temperature) T c and non-zero for T < T c as the system enters an ordered state. We will see later, 
that m oc (T c - T) 13 . In the case of a fluid, the order parameter is identified as the the difference 
between liquid and vapour density which becomes nonzero as T < T c . The analogy between mean- 
field thermal critical phenomena and percolation is not complete since the critical exponent (3 for 
the order parameter is \ for thermal phase transitions and unity for percolation, but nevertheless, 
the qualitative similarities are clear. 

In both the thermal phase transitions mentioned above and percolation, the order parameter 
goes to zero continuously as one approaches the critical point. Such phase transitions are called 
continuous phase transitions or second-order phase transitions. If instead the order parameter 
jumps to zero, one has a first order phase transition. Such transitions can occur in more complicated 
percolation theory situations, like "bootstrap percolation" where on a square lattice a site remains 
occupied only if three or four of its neighbours are still occupied [2] . 

1.5 Cluster Number Distribution 

Unfortunately, it is not possible to obtain an exact form for the cluster number n s {p) in d > 1 or the 
Bethe lattice since, unlike in Id, there are a very large number of different ways in which clusters 
can arrange themselves. Even for relatively small cluster numbers in 2d we run into difficulties. 
Exact enumeration of small cluster "animals" has been tabulated. In a 2d square lattice a computer 
has counted a total of 68,557,762,666,345,165,410,168,738 different cluster configurations with 
s = 46 [3]. Nevertheless, we are still able to account for the general scaling behaviours of certain 
quantities that characterise the clusters. 

Definition 13 The perimeter or external boundary t of a cluster is the number of empty nearest 
neighbours of the cluster. 

Using the concept of a perimeter, we are able write down the general expression for the cluster 
number n s (p) as 

n s (p) = Y,9s,tP s (l-p) t , (1.16) 
t 

where we have introduced g s j, the number of different lattice configurations with size s and perime- 
ter t. 

In Id, all clusters have two perimeter sites and thus 




1 for t = 2 
otherwise 



implying n s (p) = p s (l — p) 2 . 

Exercise 9 Consider clusters of size s = 3 in 2d on a square lattice. Verify, that g s j is only 
non-zero for t = 7 or t = 8 and that g s=3 t=7 = 4 and g s=3 t=8 = 2, implying n s= 3(p) = 4p 3 (l — 
p) 7 + 2p 3 (l -p) 8 . 
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In a Bethe lattice, however, we can show that there is a unique perimeter for a given cluster 
size s. 

Exercise 10 Show that in a Bethe lattice t = 2 + s(z — 2). 

Now we apply this formula to our general result for cluster numbers (1.16) to get: 

n s {p) = J29s,tP s (l- P y = g s , 2+s{z ^ )P s (l- P ) 2+ ^- 2) . (1.17) 
t 

For simplicity, from here on we shall work with the z = 3 case rather than the general Bethe 
lattice. For z = 3 we have p c = \. As we have said, determining g s ,2+s for large s is difficult, but 
we can avoid having to calculate g s ,2+s by considering the ratio 



n s (p) 
n s (p c ) 



(1-P) 
(1-Pc) 

(1-p) 

(1-Pc) 



2 r 



Pc 

} 

exp 



(i-p) 



sin 



1 ~Pc 

p (1-p) 



Pc (1 - Pc) 



(1.18) 



In order to proceed, define the function 

P (1 -P) _ P (1 -P) 



/(P) 



Pc(i-Pc) i(i-i) 



4p - 4p 2 



for z 



Exercise 11 Make a Taylor expansion of the function f(p) = 4p — 4p 2 around p = \ = p c and 
verify f{p) = 1 - A(p - p c ) 2 . 



Substituting the expression for f(p) given in Exercise 11, we find 



n s {p) 
n s {p c ) 



(1-P) 



(1-Pc) 



exp( 



(1.19) 



with 



In 



(l"4(p-i 



for p 



(1.20) 



In the Bethe lattice, a = 

We now have a very simple exponential decay for the ratio of cluster numbers. This is to say 
that large clusters with s » are very rare indeed. This explains why is known as the cutoff 
cluster size. In order to investigate the cluster number distribution in more detail, we have to 
consider the term n s (p c ) appearing in the denominator. 

From the analysis on the average cluster size, we know that, 



s( P ) = E 



s 2 n s (p) 
!E~i sn s {p) 



oo 



for p^ p c 



(1.21) 



However, the denominator remains finite at p = p c (why?) so the numerator must diverge for 
p = Pc, that is, 

oo 

^s 2 n s (p c ) = oo. 
s=l 
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Generally, we have that 

convergent for a < — 1 



s=l 



divergent for a > — 1. 



Thus a power-law decay 

oo oo 

n s (p c ) oc s~ T s2n s(Pc) = X! s2_T 

s=l s=l 

would imply a divergence of the average cluster size if the critical exponent r < 3. 

However, see Eq.(1.15) and problem (8), we also need p c = Y^^i sn s(Pc) oc J2^=i sl ~ T < 00 
implying that r > 2. Therefore, we arrive at a general form for the cluster number distribution in 
the Bethe lattice: 

n s (p) oc s~ r exp( ), for s > 1 (1.22) 

with 

S£(x\p c -p\~* forp^p c , (1.23) 

where a = \ and 2 < r < 3. 

Actually, we can determine r as it is not independent of the critical exponents 7 and a deter- 
mined previously. Let us evaluate the average cluster size S(p) assuming Eq.(1.22) is valid for all 
s. Neglecting the denominator which for p — ► p c will approach p c and thus be a constant we find 



s (p) K J2s 2 n s (p) 



s=l 

00 

oc y^s 2 ~ T exp( ) 



4- 



k, l s 2 T exp( ) ds 

Jl S£ 

= / (zs£) 2 ~ T exp(— z)s{. dz with z = s/s^ 

J l/s e 

/ z 2 ~ T exp(— z)dz 

Jl/s ( 
rco 

~ T I z 2 ~ T exp(—z)dz for p — > p c where st — ► 00 
Jo 

= s|" r r(3-r) (1.24) 
where T is the Gamma function 

T(x) = / z^expi-z) dz. 
Jo 

Exercise 12 Show that the Gamma function T defined by the integral above satisfies the recursion 
relation T(x + 1) = xT(x) and deduce that T(x + 1) = x\ for x = 0, 1, 2, . . . where 0! = 1. 

Using S£ oc (p — p c )~ 2 , we see that 

S{p)tx{p-p c ) 2T - & <x{p-p c y\ (1.25) 
where we have used Eq.(1.14). Thus we must require that 

2r-6 = -l^r=^. (1.26) 
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In summary, for the Bethe lattice 

n s (p) oc s _ 2exp(— — ) s S> 1 

s% oc \p-pc\~ 2 ioxp^pc. 

The above procedure is a general way of deriving a relationship between critical exponents, as we 
shall see shortly 
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Figure 1.9: (a) The cluster number distribution n s (p) = 

different values of p approaching p c = 1/(2— 1) = 1/2. The vertical lines indicate the cutoff cluster 
size s^(p) for the various p values. The cluster number distribution is well approximated by a 



power law for s <C while is decays rapidly for s 3> S£. (b) By plotting 
s/s^ s • A(p c — p) 2 , all the data collapses onto the scaling function e~ x . 

Now, we postulate a scaling ansatz (form) valid for all p and large s 
Ansatz 1 



(l-Pc) 



s T n s (p) versus 



or equivalently 



n s (p) oc s T exp(- 
S£ oc \p-p c \~° 

n s {p) oc s~ T exp(-s\p - p c \v) 



s > 1, 
P^Pc 

s>l and p ^ p c 



(1.27) 
(1.28) 



Note the role of as a cutoff cluster size, where n s (p) is characterised by 

n (p) oc / f ° r S <<C S « (1.29) 

s 1 decays rapidly for s ~> 

representing the crossover from a behaviour of "critical clusters" (power-law distributed) to that 
of non-critical clusters. Note also that n s (p c ) = s~ T as = oo at p = p c . 

We can, however, immediately identify two problems (a third problem will surface shortly) with 
the scaling ansatz 1. 
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(1) The scaling ansatz 1 fails to describe the Id case. We found previously that 



n s ip) = (1 - P) 2 exp( ) = {p c - p) 2 exp( ), 

which is not a special case of our ansatz 1 since, instead of a power of s, we have (p c — p) 2 
in front of the exponential. You might be very tempted to say that r = 0, but that does 
actually not solve the problem, as we shall see below. 



(2) For fixed s 



(p)=T,9s,tP s (l-p) t (1.30) 



n _ 

t=i 



is a polynomial in p with a finite number of terms. Thus all derivatives of n s (p) with respect 
to p will remain finite for all p. However, using the scaling ansatz 1 for p > p c n s {p) = 
As~ T exp(— s(p — p c )~) we find, keeping s constant 

^ = As -r eM _ s{p _ p c )*)I(p _ p c )*-l (1.31) 

dp a 
and by differentiating once more 



d 2 n s {p) _ ^ 
dp 2 



vxp(-s(i> - p c )°) [\p-Pc) 1 " ^) + exp(-s(p - p c )*)^ - 1 ) ip - /v) 



In 2d, e.g., a = 0.4 <^ ^ — 2 = 0.5, so both the first and second derivative of n s (p) w.r.t. p 
remain finite as p — > p c . The third derivative of n s (p), however, will contain a term with the 
factor (p — p c ) ~ 3 which will diverge as p — ► p c contrary to the derivatives of Eq.(1.30). 

We shall later resolve these two problems, so there is no need for you to worry too much at the 
moment. 

Let us return to the question of deriving scaling relations assuming the scaling ansatz above for 
the cluster number distribution. Following the previous calculation of the average cluster size S(p) 
and noting that the Gamma function is just a number, we find 

S(p) cc s|~ T oc \p - p c \^~, ioip^pc, (1.32) 

S£<x\p-p c \~v ioip^p c . (1.33) 
S{p) oc \p -p c \~~ f ioxp^pc, (1-34) 



as 



However, by definition 
implying the scaling relation 



7=—- (1-35) 
a 

We notice again that to ensure the divergence of S(p) for p — > p c we must have r < 3 in order to 
leave 7 > 0. Also, we recover problem (1) mentioned above because in Id percolation theory 7 = 1 
and <T = 1 implying r = 2. However, there is no s~ T term (yet!) in our expression for the Id cluster 
number distribution n s (p) and obviously the suggested "solution" of assuming r = is not valid 
either. 
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Now we focus our attention on the percolation strength P(p), starting from Eq.(1.15): 



P(p) = p~Y^ sn s (p) 
s=l 

oo 

= Pc-^sn s (p) + (p-p c ) 
s=l 

oo oo 

= ^2sn s (p c ) - ^2sn s (p) +0(p-p c ) (0(p - p c ) = term of order (p - p c )) 



s=l s=l 

oo oo 



OC 



OC 



s 1 T — s 1 T exp( ) + 0(p — p c ) using the scaling ansatz 1 

S=l 5=1 S Z 

E^ T {l-exp(--))+0(p- Pc ). 



s=l 

oo 



S 
S 



E» 1_T l-exp( ) As, (1.36) 



s=l 



where we have dropped the term of order (p — p c ) and introduced As = 1. When p — ► p c , the cutoff 
cluster size s^ — > oo. Thus the term {1 — exp(— — )} ~ for small s and the main contribution to the 

integral will come from large s of the order of s^. To solve (1.36) we need to replace the summation 
by an integral. Going from a discrete sum to a continuous sum does not, of course, create a 
mathematical identity, unless one takes certain limits. However, when the main contribution from 
the sum is for large s, it is legitimate to approximate the sum by an integral ^ As ~ / ^ s > as we 
are only interested in the the general scaling behaviour of the system. Thus 

P(p) oc J™ s l ~ T jl -exp(--^)j ds. (1.37) 

One might proceed by integration by parts. You will then discover that only one of the terms are 
non-zero. However, we will take another approach. We differentiate Eq.(1.37) w.r.t. 1/s^ to find 

— — r— oc / s exp( las 

d(i) Ji PV V 

oc s? _T / z 2 ~ T exp(— z) dz substituting z = s/s%. 

Jl/s ( 

As we are ultimately interested in the limit p — > p c where the cutoff cluster size diverges, we replace 
the lower limit 1 /s^ with zero and thus 

P{p) oc (^j (1.38) 

so integrating we find 

P{p) = Cl + C2 (j-^j (1.39) 
At p = pc, s^ = oo. But, P(p c ) = implies that C\ = 0. So, 

Pip) oc C2 (l) , 

oc \V-Vcf^- (1-40) 
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Hence, we have derived a new scaling relation 



By definition 



P=—. (1-41) 



P(P) = | ? y, P ~ \ (1-42) 



Now we find a third problem with our ansatz 1. 

(3) In the derivation above we did not assume p > p c . Thus the ansatz 1 also predicts an infinite 
cluster for p < p c , which is of course not correct. 

Reflecting upon the three problems, we realise that they can all be resolved, as we shall see 
shortly, by avoiding the argument s\p — p c \ * in the scaling ansatz. Thus let us assume the following 
general scaling form 

Ansatz 2 

n s (p) = qos~ T f[qi(p - Pc)s a ] s > 1 and p -> p c . (1.43) 



The scaling ansatz 2 has the following properties: 

• go and q\ are proportionality factors and p c the critical occupation probability. They all 
depend on the lattice details. They are non-universal. 

• t, a are the critical scaling indices which are independent of p,p c and the lattice structure. 
They are universal. 

• /, the scaling function, is independent of p,p c and detailed lattice structure. It is univer- 
sal. The precise form of the scaling function / = f{z) has to be determined by computer 
simulations or other numerical methods since the function / is not predicted by the scaling 
assumption. Typically, however 



i 



f i \ — J constant \z\ <C 1 44> \p — p c \s a <C 1 O s <C \p — p c \ ° oc s^ ,-. . ... 

~~ \ decays fast j^j > 1 4^ s > S£. 

that is, the scaling function f(z) reaches a constant value (except in Id percolation theory) 
for \z\ <C 1 <^ s <C S£, and decays rather fast for \z\ » 1 44> s » s^. Thus, the role of s^ as a 
cutoff and crossover size is maintained. Later, we shall look more carefully at the crossover 
properties. 

• r, a and / depend on dimensionality. 

• By plotting j^s T n s {p) versus (p—p c )qis a , all the data will fall on the universal curve outlining 
the scaling function /, see Fig. 1.9. This phenomenon is known as data collapse. 

Problems 7 and 8 highlights some of the points above. 

Let us check whether the general scaling ansatz 2 really resolves the three problems associated 
with the former scaling ansatz 1. 
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(1) The Id cluster number distribution 

n s (p) = {Pc -f>) 2 exp(-— ) 

~ {Pc - p) 2 exp(-s(p c -p)) ioxp^pc 
= s' 2 (s(p c -p)) 2 exp(-s(p c -p)) 
= s- 2 f[s°{p c -p)} 

(1.45) 

which is of the form of the general scaling ansatz 2 with r = 2 and a = 1 and f(x) = 
x 2 exp(— x). Thus the scaling relation Eq.(1.35) is indeed fulfilled. Note, however, that unlike 
in higher dimensions, the scaling function increases like x 2 for x <C 1 but that it indeed decays 
rapidly for x S> 1, see Fig. (1.5). 



(2) Keeping s constant we find 

qos- T f (k) [<n(p-PcK}qiS 



d k n s (p) _ T {k) _ s ff]n k a ak 



Qpk 

which is not divergent for p — > p c assuming that / is analytic, that is, all derivatives of f(z) 
with respect to z are finite everywhere, also at z = 0. 

(3) Let us calculate the strength of the infinite cluster. 

oo 

-Pi?) = ^Z sn s{p)~P 

8=1 

OO 

= SU s(p) ~Pc + (Pc - P) 

3=1 

OO 

= Y} sn s{p) - sn s (p c )} +0(p- p c ) using £ sn s (p c ) = p c 

s=l 

/oo 
qos 1 ^ if[qiip-Pc)s a ]-f[0])ds 

Now substitute z = q\{p — p c )s a . Notice that if p > p c the argument z is positive while 
for p < p c the argument is negative. The cluster size s, however, is always positive so 
s= \z\°(qi\p — Pc\)~~ implying that ds = ^\z\~~ 1 (q\\p — p c \)~°dz with dz positive. Thus we 
find 

-Pip) = Jqo\z\ 1 -^iqi\p-p c \) T -^[fiz)-f(0)}^\z\^- 1 (q 1 \p-p c \)-^dz 

= qoiqilp-Pcl) 1 ^ j\z\ 2 -^- 1 [fiz)-f(0)]^dz. (1.46) 

Now it is important to realise that the limits of the integral will depend on the value of p 
since z = q\{p — p c )s a . If p > p c , the upper limit will be oo (lower limit p — p c approaches as 
p — > Pc) while for p < p c , the lower limit will be — oo (upper limit approaching as p — > p~). 
Thus if 

1 ^ _ I constant / for p > p, 
1 d I for p < p c 



I 



z \^-\f{ z ) - f(o)}- dz = { ;; U1 ™ ^ u (L47) 



we have solved the third problem as well, implying 

P(p) oc 



ip — Pc) 13 for p — > p+ 
for p < p c 
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with 

a 

Exercise 13 Discuss the qualitative behaviour of the scaling function f(z) as a function of z = 
Qiip ~ Pc)s a in order to satisfy the constraint given in Eq.(l.J^l). 

Problem 8 will ask you to calculate the average cluster size S(p) assuming the general scaling 
ansatz 2. 

1.6 Cluster Structure 

Till now we have only considered the distribution of cluster sizes and have come up with some useful 
scaling laws and scaling relations. Now we turn our attention to the geometry of the clusters. We 
shall do so by mainly studying the fractal geometry, which contains information about the density 
of the clusters on different length scales. 

1.6.1 Cluster Radius and Fractal Dimension 

Let Yi denote the position of the ith occupied site in a cluster of size s. 
Definition 14 The centre of mass of an s-cluster r cm = ^ J2t=i r i- 
Definition 15 The radius of gyration R s of a given s-cluster is defined by 
R 2 S = average square distance to the centre of mass 




We want to prove that 

R 2 S = i average square distance between two cluster sites = l r « — r j\ 2 - (1-48) 

By definition, the lefthand side of Eq.(1.48), 
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Let us rewrite the righthand side of Eq.(1.48) 



R ' ~ 2^ 



s s 



i=l j=l 



1 1 

2^ 



Y s Y i ■ Y i + Y ■ 

i=l j=l 




1 S 1 

S i=l S 
1 S 

1=1 

The radius of gyration is a more useful quantity to work with since in many situations (like in 
Polymer Science) we have to deal with more complicated structures than straight lines, squares or 
spheres. 

We also remember that the correlation function ^(r) is defined as the probability that a site at 
a position r from an occupied site belongs to the same finite cluster. We have the identity 

EsW = %) (1-49) 

r 

because the sum is the average number of sites to which an occupied site is connected, see Problem 
2. 

Definition 16 The correlation length £ is defined as 

* " E r <7W • (L50) 

The correlation length represents some average distance of two sites belonging to the same cluster. 

For a given s cluster, 2R 2 S is the average squared distance between two cluster sites. The 
probability of a site being part of an s cluster is sn s (p), and it will be connected to s sites (if we 
include self-connection). The corresponding average over 2R 2 S gives the squared correlation length 

^ " J2 s s*n s ( P ) ■ (L51) 

The correlation length £ is the upper cutoff of the radius of those clusters which contribute to 
the mean cluster size near the percolation threshold. It is expected that £ diverges as p — > p c , like 

Z<x\p-Pc\-"- (1-52) 

We have introduced a new critical exponent, u, and we wish to determine how it is related to r 
and a. For 2d percolation, plausible but not rigorous arguments give v = |, in excellent agreement 
with numerical results. In 3d, v is somewhat smaller than 0.9, whilst for Bethe lattices one has 
v = ry, analogous to the mean-field theories for thermal phase transitions. 

As we have already seen, many quantities diverge at the percolation threshold. Most of these 
involve sums over all cluster sizes s; their main contribution comes from s of the order s^ oc 
\p — p c \~~ ■ Now we see that the correlation length £, which is also one of these quantities, is simply 
the radius of those clusters which contribute the most to divergence. This effect is the foundation of 
scaling theory. There is one and only one length £ dominating the critical behaviour in an infinite 
lattice. 
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1.6.2 Finite Boxing of Percolating Clusters 



We now want to find out how the R s varies with s at the percolation threshold p = p c where £ = oo. 
Let M{L) denote the mass of the percolating cluster within linear distance L. If the percolating 
cluster was a compact object, then it would be of the form M(L) oc L 2 . However, at p = p c , 
the percolating cluster is a fractal object. Figure 1.10 displays the size of the largest cluster 




Figure 1.10: The size of the largest cluster 5^ at p = p c as a function of the lattice size L. 

(the percolating cluster) as a function of lattice size L in a double-logarithmic plot. The data are 
consistent with the existence of a fractal dimension D such that 



M(L) oc L D . 

For the 2d case D = §§ < 2 and in Id, D = 1. 

Having obtained the relationship in (1.53) it is natural to assume 

s oc for s 3> 1 at p = p c . 

Thus, 



(1.53) 



R s oc so 

Inserting Eq.(1.55) into Eq.(1.51) we have 



for p -> p c . 



(1.54) 
(1.55) 



e 2 oc 



oc 



J2 s ^ 2+ on s (p) 
J2 s s 2 n s (p) 



\p-p ( 



IT-3- 



.121 



where we have used 



ip-p c r _l_ ~ 
« ip-p c r 2i/ , 



^2s k n s (p) (x\p-p c 

8=1 



(1.56) 
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Therefore we end up with the scaling relation 



Da 



44> D 



1 

va 



(1.57) 



Note that the behaviour at p = p c and the behaviour for s <C s^, or R s <C £ are indistinguishable. 
Therefore, £ is the crossover length between a critical and noncritical behaviour. On length scales 
much less than £ we have n s oc s~ T and R s oc so but for length scales much larger than £ the 
scaleless behaviour disappears. Note also, that 



S£ oc £ D oc |p — Pc| 1/0 oc |p- 



confirming once again the scaling relation above between D, v and a. 

If, however, p 3> p c , the fractal dimension D is expected to be equal to the Eucledian dimension 
d of the lattice. So, 

D = d for p > p c . (1.58) 

Clusters above p c are not fractals but 'normal' objects with D = d, provided s > or equivalently 
R s » £. Thus or £ sets the cluster size and linear scale where there is a crossover from fractal 
to non-fractal behaviour. 



Exponent 


Id 


2d 


3d 


4d 


5d 


6d 


Bethe 


a 


1 


-2/3 


-0.62 


-0.72 


-0.86 


-1 


-1 


P 





5/36 


0.41 


0.64 


0.84 


1 


1 


7 


1 


43/18 


1.80 


1.44 


1.18 


1 


1 


V 


1 


4/3 


0.88 


0.68 


0.57 


1/2 


1/2 


a 


1 


36/91 


0.45 


0.48 


0.49 


1/2 


1/2 


T 


2 


187/91 


2.18 


2.31 


2.41 


5/2 


5/2 


D{p = p c ) 


1 


91/48 


2.53 


3.06 


3.54 


4 


4 



Table 1.2: The critical exponents for the percolation theory problems in dimensions d = 1, 2, 3, 4, 5, 6 
and in the Bethe lattice. 



1.6.3 Mass of the Percolating Cluster 

Now we wish to see how the scaling relation of the mass M{L) changes for the two different cases. 

(1) For L <C £: (always the case at p = p c , where £ = oo) The percolating cluster appears fractal, 

implying M(L) oc L D . 

(2) For L » £: The cluster appear homogeneous 

M(L) = (No. of lattice sites) x (Prob. site belongs to percolating cluster) 
= L d P(p) 
= L\p-p c f 

since £ oc \p — p c \ v - (1.59) 



Now consider L k, £, that is, we match the observations above by substituting L with £: 

M(L) oc L D oc £ D , (1.60) 
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and 

M(L) oc L d ^^ oc £ d ~v. (1.61) 

Therefore, we can state that 

D = d-P, (1.62) 

which is know as a hyperscaling relation because the Euclidean dimension d enters in the scaling 
relation. 

Since the percolating cluster has a constant density for L 3> £, it is natural to divide the system 

into boxes of linear size £. In cZ dimensions, the total volume L d will be divided into (j^j boxes. 

Since the cluster inside each of these boxes of size £ d has a mass of order £ D , the total mass of the 
cluster is given by 

M(L,0=(^)\ D = e- d L d , (1.63) 

which is, of course, equivalent to P(p)L d . 
In summary, we have 

I L D L«f 



which can be written in an alternative form of 

= T, D mi 



M(L,0 = L D m&, (1.65) 



where the scaling function 



, N f constant x = h <C 1 , 
™(*H^ x=|»l. (1-66) 



The argument x of the scaling function is a dimensionless number, namely the ratio of the two 
length scales j. No other length scales play a role. The scaling function describes a crossover 
from fractal behaviour at length scales L much smaller than the correlation length £ to uniform 
behaviour at length scales L much larger than the correlation length £. In Nature, there are many 
examples of crossover phenomena from one type of behaviour to another. Take as an example the 
table in front of you. On large length scales, the surface of the table is smooth, but going to very 
small length scales, the surface will become very rough and probably fractal. The mountain range 
is yet another example. Far away, i.e., on large length scales, it looks pretty smooth, but, as you 
are well aware, if you are in the middle of the mountain range, on small length scales, it is pretty 
rough and rugged. 

Exercise 14 Give one or two more examples of such crossover phenomena in Nature. 
The density has the form 

P( L ,0= Ld ^D-d L> £ = L P(j), (1-67) 
with a scaling function 

, f constant x = h <C 1 

**> - { ^ x = | » i. < L68 > 

To be aware that the density actually decreases with the length scale when L < £ can be of 
immense practical importance. Imagine that you have just been hired by an oil company and you 
were given a sample of dimensions 0.1 x 0.1 x 0.1 m 3 of the porous medium where the oil resides and 
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Figure 1.11: The density of sites belonging to the percolating cluster within a region of linear size 
L. In region I, the ratio L/^ <C 1 and in region II, L/£ 3> 1. Region III is related to the fact, that 
the lattices used are finite. 

had to estimate how much oil the company could expect to recover in the oil field of dimension, 
say, 100 x 100 x 100 km 3 . The percolation model is used as a model for the distribution of oil 
and gas inside porous rock in oil reservoirs. The empty sites model hard rock while the occupied 
sites are pores filled with oil. The average concentration of oil is p. Now, if p < p c there are 
only finite clusters, so let us assume p > p c . The strength P(p) will be the probability of drilling 
into a percolating cluster, but how much oil can you recover? You would head for an immediate 
dismissal if you argued as follows: Measure the density of oil in the sample p{L\ = 0.1m). The 
mass of oil to be recovered is M = p(L± = 0.1m) • (100000m) 3 . The argument is wrong simply 
because the density of oil at length scales L 2 = 100 km is not given by the density of oil at length 
scales L\ = 0.1 m, see Fig. (1.11). The proper way to argue would be the following: measure the 
density of oil in the sample p(L\ = 0.1m) = CL®~ d . The density of oil at the reservoir length scale 
p{L 2 = 100km) = CL 2 D ~ d with the same constant C. Taking the ratio we find 



since in 3d D — d « 2.53 — 3 = —0.47, see table 1.2. The former estimate is wrong by a factor 1000. 

1.7 Fractals 

The percolating cluster at p = p c is an example of a random fractal. When p / p c large clusters 
appear fractal on length scales up to the correlation length £ oc \p — p c \~ u , roughly speaking. In 
mathematics, we can also encounter deterministic fractals as for example the Sierpinski carpet and 
the Sierpinski gasket. 

For a geometrical fractal, the mass scales with the linear size raised to the fractal dimension 



p(L 2 = 100km) 




« 1.5 • 10~ 3 p(Li = 0.1m), 



M(L) = L u ^D = 



log M(L) 
logL 
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Figure 1.12: Two examples of deterministic fractals, (a) The Sierpinski carpet. The first three levels 
in the algorithm for constructing the Sierpinski carpet are shown. The linear scale is enlarged by a 
factor 3 with the condition that an occupied square (white) is replaced with with 3x3 squares in 
which the centre square is empty (grey) while an empty square (grey) is replaced with 3x3 empty 
squares, (b) The Sierpinski gasket. A similar algorithm generates the Sierpinski gasket. 



Let n denote the iteration number. By inspection of the Sierpinski carpet, M = 8 n and L = 3 n , 
implying 

log8 n log 8 
D = -2 — = -°- « 1.893. 
log 3™ log 3 

In the case of the Sierpinski gasket, M = 3 n and L = 2™, yielding 

D=pl* 1.585. 
log 2 

Note that if the length scale is rescaled by a factor b, then 



b' \bj " b D ~ b D 

so that the mass is reduced by a factor b D . Using this observation as the characteristic of fractal 
behaviour for random fractals we can easily find, say, the scaling of the mass of the percolating 
cluster M(L,£) as a function of L and £. Let us rescale the length scales by a factor b, that is, 
L — > j and £ —>■ | (no other length scales are present). Thus 

M(L,0 = b°M(j,i) 

= b°b°M^,l 2 ) 



= b Dl M(- 
[ b l, b l) ' 



(1.69) 
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If we are at the critical point p = p c , £ oc \p — p c \ v = oo, but then ^ = oo V/. If we stop the above 
renormalisation procedure when b l = L we find 



M(L,£ = oo) = (V) D M(^,i = oo) 



= L D M(l,oo) 
= L D . 

(1.70) 

This proves that self-similarity is mathematically expressed by power-law behaviour. Now, assume 
p > p c with £ oc \p — p c \~ v < oo with 

We now have to consider two cases depending on the ratio L/£. 

(a) If L/£ <C 1 44> L <C £ (which is airways the case at p = p c ) we stop the iteration process when 

we reach the smallest length scale b l = L and we "recover" 

M(L, = L D M(1, |) = L D M(1, oo) oc L D . 

(b) If L/£ » 1 L » £, we stop the iteration process when b l = £ and we find 

= e D M(|, i) oc e D (|) d = 

since for a uniform system M(j, 1) oc . 
Exercise 15 Discuss why the system looks uniform when £/b l = 1. 

1.8 Finite-size scaling 

Let us consider the percolation problem for p ^ p c where the correlation length is finite 

£ oc \p-p c \~ u =>• \p-p c \ oc £ - -. 

The strength of the infinite cluster for p > p c can then be expressed in terms of the finite correlation 
length 

P(p) OC(p- Pc) 13 OC 

and similarly for the average cluster size 

S(p) oc \p-pc\~~ 1 oc 

These are the results in infinite lattices with L = oo. Now, what happens at finite lattice sizes 
when L < oo? Well, nothing happens as long as L > (. The only relevant length scale is the 



29 



correlation length £. However, when L <C £, the length scale will be set by L. This we have already 
encountered before, we just didn't stress that point. Better late that never: 



where the scaling function 



P(L,£) = p(L,0 
M{L,i) 



/(>) 



r«/(f) 



4) = ((*)"" H 1 



^ 1 constant ^ S> 1. 




Figure 1.13: (a) The strength P{p) measured in a 2d square lattice for finite system sizes as a 
function of occupation probability p. When p is far away from p c , L £ and the strength is 
determined by ■ However, with p closer to p c , the finite-size effect is seen and P(p) decays with 

system size as L~~ . (b) The average cluster size S(p) measured in a 2d square lattice for finite 
system sizes L as a function of occupation probability p. When p is far away from p c , L>^ and 

7 

the average cluster size is determined by the correlation length £ *> . However, with p closer to p c , 
the finite-size effect is seen and S(p) increases with system size as L~ . 



Similarly for the average cluster size, we would expect 

S(L,£)<x 



Generally, if we have a quantity 

X oc \p — p c \~ x °c £~ for L 3> £ 
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then we would expect 



X(L,0 oc 



i Li L«£ 

l » e 

(1.71) 



or, in terms of occupation probabilities 

X(L,p) = \p - p c \~ x x 2 (Lv (p - Pc )). 



The important message we get is that by studying finite-size scaling, that is, studying X as a 
function of (finite) system sizes L at p = p c , we can extract the exponent 



1.9 Real space renormalisation in percolation theory 

Imagine a percolation lattice with n sites. The total number of different microscopic configurations 
is 2 n as a site can be in one of two states - occupied with probability p or empty with probability 
1 — p. Table (1.3) shows the number of microscopic configurations in a two dimensional square 
lattice of linear size L. 
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16 
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6 • 10 i4 
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2 • 10 24 


1Q 3U 


10« v 
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Table 1.3: The number of microscopic configurations 2 n as a function the linear size L of a two- 
dimensional square lattice where n = L 2 , the number of sites in the lattice. For lattice size L = 17, 
the number exceeds the estimated number of atoms in the universe « 10 81 . 

In principle, if we wanted to calculate the average of a quantity {A) (say the average cluster 
size), we should evaluate the quantity Ai associated with each configuration i and then weight the 
quantity Ai with the probability pi (which would be of the form p m (l — p)( ra_m ), where m is the 
number of occupied sites in that particulate configuration) of the system to be in that microscopic 
state, that is, 

<a> = ^p^ = 5>*. (1.72) 

This procedure, however, soon becomes an impossible task due to the numbers of different 
states involved. With today's computing power one might do, say, a 7 x 7 system, but an 8 x 8 
system would take approximately 40 years on a tera flop (10 12 calculations per second) while the 
time required to do a 10 x 10 system would exceed the age of the universe which is approximately 
15 • 10 9 years. 

Real space renormalisation is a coarse graining procedure where we systematically reduce the 
number of degrees of freedom (sites) in order to break a large problem down into a sequence of 
smaller and more manageable stages by eliminating fluctuations on scales less than a given length 
scale b. The effect will be to reveal the large scale behaviour of the system. 

Before we introduce the real space renormalisation method, let us make some observations 
related to the phase transition in a percolation model. 

(a) The empty state at occupation probability p = is self-similar with an associated correlation 
length £ = 0. 
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(b) The fully occupied state at occupation probability p = 1 is self-similar. The associated corre- 

lation length £ = 0. The ordered state has an order parameter P(p = 1) = 1. 

(c) The states at the critical occupation probability p = p c are also self-similar and the correlation 

length £ = oo. The infinite cluster is fractal and looks alike on all length scales and we have 
clusters of all sizes from s = 1 to s = oo. 

Loosely speaking, the correlation length £ sets the scale of the largest fluctuation from the 
"averaged" state. For p / p c , the correlation length is identified with the linear size of the cutoff 
cluster size s%. 

When p < p c the upper cutoff of the linear size of the fluctuations away from the empty state 
is given by £ and when p approaches zero, £ approaches zero. At p = 0, there are no clusters left 
at all and £ = 0. 

When p > Pc, the correlation length can also be identified with the upper cutoff in the linear 
size of the largest holes in the percolating (spanning) cluster which indicates how far away you are 
from the fully occupied state. When p approaches one, £ approaches zero. At p = 1, there are no 
holes left and £ = 0. 

At the phase transition point p = p c , the correlation length is infinite. Thus there are fluctua- 
tions from the smallest length scale up to infinity. 

The real space renormalisation technique is based on the so-called block site (spin) technique in- 
troduced by Leo Kadanoff and later formalised into the renormalisation group method by Kenneth 
Wilson, who received the Nobel price in 1982 for "his theory for critical phenomena in connection 
with phase transitions", see, for example, the website on physics laureates http:/ /www. nobel.se/physics/laureates/ 

The renormalisation group method has three basic steps: 

1. Divide the lattice into blocks of linear size b (in terms of the lattice constant) with each block 
containing a few sites (spin). 

2. Next, the coarse graining procedure takes place. The sites in the blocks are averaged in some 
way (to be specified more precisely shortly) and the entire block is replaced by a single su- 
per site (spin) which is occupied with a probability according to the renormalisation group 
transformation p' = R b {p). 

In the combined procedure 1 and 2, one should keep the symmetry of the original lattice such 
that we can repeat the coarse graining procedure again. The result of these two operations are 
to create a new lattice whose fundamental spacing is b times as large as the original lattice. 

3. Restore original lattice constant by rescaling the length scales by the factor b. 

These 3 steps define a renormalisation group transformation R b , where the number of degrees 
of freedom is reduced drastically. Note that it is called a group because of the property that when 
applying the transformation to a configuration 

R b2 (R bl ({ Si })) = R blb2 ({ Si }) (1.73) 

but it is of course not a group in a strictly mathematical sense because no inverse transformation 
exists since we are reducing the number of degree of freedom. 

The effect of the coarse graining procedure in step 2 is to eliminate from the system all fluctu- 
ations whose scale is smaller than the block size b. Any small scale fluctuations of the sites over 
a range of less than b lattice units will be smeared out. This is somewhat similar to the process 
taking place if you view the system through a lens out of focus - all the smaller features are blurred 
but the large scale features are unaffected. 

Let us clarify the 3 steps by looking at a simple example. Consider a 2d square lattice with 
9x9 sites, that is, 81 degrees of freedom and a total of 2 81 = 2 • 10 24 different configurations. 
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Figure 1.14: The three renormalisation group steps illustrated in a 2d square lattice with b 
The majority rule defines the renormalisation group transformation. 



3. 



1. Divide the square lattice up into small 3x3 blocks with a total of 9 sites (6 = 3.) 

2. We now have to specify the course graining procedure which will define the renormalisation 
group transformation. One possibility would be to use the so-called majority rule which 
states that if a majority of the sites in the block is occupied, then the super site should be 
occupied, otherwise empty. The probability, that a majority of sites are occupied in the 3x3 
block then defines the renormalisation group transformation 

jf = R b (p) =p 9 + K 9AP 8 (1 -p) + K 9 , 2 p 7 (l - pf + K 9:3P 6 (1 - pf + K 9A p\l - p)\ 

where Kg^ is the number of different ways to place k empty sites into the cells of 9 sites. 
The super site should be occupied with probability p' = Rb(p) and empty with probability 
l-p> = 1-R b (p). 

3. Rescale the lattice of super sites by a factor 3 to restore original lattice constant. 
Let £ denote the correlation length in the original lattice. Assume we are close to p c . Then 

£ = constant \p - p c \~ v (1.74) 

We denote the correlation length in the new renormalised lattice by Since the length scales 
have been rescaled by a factor b (step 3) 

C = I (1-75) 



constant \Rb(p) —p 



_ v constant \p — p c \ 

cl 



where Rb{p) is the renormalisation group transformation that determines the new occupation prob- 
ability in the rescaled lattice associated with the new (smaller as b > 1) correlation length. 
Simple rearranging implies that 



f\R b ( P )-Pc\Y u 1 



so the critical exponent v is given by 



V \p~Pc\ J b 
log b 



(1.76) 



(1.77) 



The basic idea is that at the critical point p c we have self-similarity, that is, 

= e (us) 
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Figure 1.15: When renormalising the length scales by a factor b, the correlation length decreases 
unless it was or oo initially. Associated with a decrease in the correlation length is a flow in 
occupation probability space either towards p = 0, if initially p < p c or towards p = 1, if initially 
p > p c . The fixed points of the renormalisation group transformation are associated with the trivial 
self-similar states at p = or p = 1 or the nontrivial self-similar states at p = p c . 




This equation can only be consistent with Eq.(1.76) if 

?={ °° (1-79) 

This is intuitively clear as £ = is associated with either p = (empty lattice) or p = 1 (fully occu- 
pied lattice) the two trivially self-similar cases and £ = oo with the critical point p c . Furthermore, 
we can also conclude that 

£ = constant \p — p c \~ v = £' = constant \Rb{p) — pd\~ v (1.80) 

which implies that the phase transition £ = oo is identified with (one of) the fixed points p* of the 
renormalisation group transformation, that is, the solutions of the equation 



Rb(p*) = P*. (1.81) 

Thus 

\Rb{p)-Pc\ \R b { P ) - R b { Pc )\ 



\p - p c \ \p - p, 

dRb(p c 



dp 

leaving, see Eq.(1.77) 



for p -+ p c (1.82) 



1.9.1 Renormalisation group transformation in Id. 

Consider a renormalisation group transformation in Id. Take a l<i lattice where each site is occupied 
with probability p. Divide the lattice into blocks with b sites. Let the renormalisation group 



34 



transformation be determined by the probability of having a spanning cluster. The probability of 
having a spanning cluster in a block of b sites is 

R b ( P ) = p b . 

The fixed point equation can easily be solved 

Rbivl = tf h = P *e>p* = 

If we start with an empty lattice p = 0, the renormalised lattice will also be empty, associated with 
the fixed point p* = 0. If we start with a fully occupied lattice p = 1, the renormalised lattice 
will also be fully occupied since all the blocks contain a percolating cluster, associated with the 
fixed point p* = 1. However, if we start with a lattice containing some empty sites p < 1, the 
renormalised lattice will contain even more empty sites because p' = Rb(p) = p b < p. Repeating 
the renormalisation procedure will gradually bring the renormalised occupation probability towards 
zero. Clearly, p* = is associated with £ = 0, while p* = 1 is associated with £ = oo and must 
thus be identified with the critical occupation probability p c . 

Furthermore, in order to calculate the correlation length exponent v we need to take the deriva- 
tive of the renormalisation group transformation evaluated at p* = 1 . 

—^\ P *=i=bp | p *=i = 6, 

implying 

= logfe = 1 

i0 & dp \p =1 

so in Id the renormalisation group transformation is exact. 




1.9.2 Renormalisation group transformation on 2c? triangular lattice. 

Divide the lattice into triangular cells containing three sites each. Let the probability for having a 
spanning cluster define the renormalisation group transformation. Thus 

R b (p) = p 3 + 3p 2 (l -p) = 3p 2 - 2p 3 . (1.84) 

The fixed point equation 

R b (p*) = 3p* 2 - 2p* 3 =p* <^p* = 




The two trivial fixed point is associated with the self-similar states of an empty lattice and fully 
occupied lattice. The unstable fixed point p* = 1/2 will be associated with the nontrivial self-similar 
states at p c . The critical exponent v is 

log b log a/3 log a/3 
v = TFn~\ = ; 7~„ = ; tttz = 1.355. 



l og ^M| pJ[=1/2 log(6p-6p 2 )| p * =1/2 log 3/2 



dp 

The exact values are p c = 1/2 and v = 4/3 so the renormalisation group transformation does a 
good job. 
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Figure 1.16: Renormalisation group method on a 2d square lattice according to the majority rule 
for 3x3 blocks. Shown are part of the original lattice L = 729 and the renormalised lattices 
L = 243 (part of), L = 81,27, and 3. (a) Initial occupation probability p < p* in the original 
lattice and the flow is towards the empty lattice, (b) Initial occupation probability p > p* in the 
original lattice and the flow is towards the fully occupied lattice, (c) Initial occupation probability 
p = p* in the original lattice. This is a fixed point for the renormalisation group transformation 
and consequently there is no flow. It looks like itself on all length scales. 



1.9.3 Renormalisation group transformation on 2d square lattice of bond per- 
colation. 

Consider the block defined in Fig. (1.17) which, after the coarse graining procedure, is replaced by 
two superbonds. The renormalisation group transformation is 

p' = R b {p) 

= probability to have a spanning cluster in horizontal direction 

= p 5 + p 4 (l -p)+ 4p 4 (l - p) + 2p 3 (l - p) 2 + 2p 3 (l - pf + 4p 3 (l - pf + 2p 2 (l - pf 

= 2p 5 - 5p 4 + 2p 3 + 2p 2 (1.85) 



with the fixed point equation 

R b (p*) = 2p* 5 - 5p* 4 + 2p* 3 + 2p* 2 =p* ^p* = { 




1 

1/2. 
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Figure 1.17: Renormalisation group method of 2d square lattice with bond percolation according to 
the spanning cluster rule, (a) and (b) After the coarse graining procedure, the block is replaced with 
two super bonds AC and AG. (c) The various configurations with a horizontal spanning cluster. 
The two dangling bonds DG and EH need not be considered as they do not affect the probability 
of having a spanning cluster from left to right. 



Again, the two trivial fixed points are associated with the self-similar states of an empty lattice 
and fully occupied lattice. The unstable fixed point p* = 1/2 will be associated with the nontrivial 
self-similar states at p c . For bond percolation in 2d, p c = 1/2, see Table 1.1. The critical exponent 
v is 

„= J° ( g , b = 1 1 ° g2 / = 1.428, 

which should be compared to the analytical result of v = 4/3 in 2d. Thus the renormalisation 
group transformation defined above for bond percoaltion in 2d results in an exact prediction of p c 
and a good estimate of v. 



1.9.4 Why is the renormalisation group transformation not exact? 

There are two significant sources of error which is the reason for not obtaining the exact results 

(1) Two sites that were connected in the original lattice can be disconnected and vica versa. 

(2) The renormalised system is not really a true percolation system because the probabilities of 

having bonds between supersites are no longer independent. As we shall see later, this is 
the equivalent of introducing next-nearest-neighbour couplings in the renormalisation group 
procedure in the 2d Ising model. Thus we should really introduce new parameters but we 
choose here to truncate the problem to only one parameter, the occupation probability p. 
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Glossary 



Taken from The New Physics ed. Paul Davies, Cambridge University Press. 

Coarse-graining An operation implementing some form of spatial averaging which smoothes out 
relatively small length scale configurational structure while preserving the larger length-scale struc- 
tures. 

Correlation length The correlation length £ gives a measure of the typical length scale over which 
fluctuations of one microscopic variable are correlated with the fluctuations of another. In percola- 
tion theory, it is the typical cluster diameter of the clusters which give the main contribution to 
the divergence of the second (and higher) moments of the cluster distribution. Close to a critical 
point £ oc \p — Pc| ^ for p — > p c . 

Critical phenomena The phenomena which occur in the neighbourhood of a continuous (second 
order) phase transition, characterised by very long correlation lengths. 

Critical exponents (or indices) Near a ctitical point, physical quantities are often proportional 
to a power of another quantity, such as the distance from the critical point \p — p c \. The power that 
occurs known as a critical exponent (index). 

Critical point A point in a phase diagram, where the correlation length associated with the physical 
system is, in principle, infinite. 

Fractal geometry Generalisation of Euclidean geometry suitable for describing irregular and frag- 
mented patterns such as the percolating cluster at p = p c . A noninteger fractal dimension D can 
frequently be associated with such patterns. 

Order parameter A variable such as the strength of the infinite cluster in percolation theory (or 
the magnetisation in an Ising model) used to describe the degree of order in a phase above (below) its 
critical point. In a continuous phase transition (second order phase transition), the order parameter 
goes continuously to zero as the critical point is approached from above (below). 
Phase transition A change of state such as occurs in the boiling or freezing of a liquid, or in 
the change between ferromagnetic and paramagnetic states of a magnetic solid. An abrupt change, 
characterised by a jump in an order parameter is known as first order; a change in which the 
order parameter evolves smoothly to or from zero is called continuous or second order. 
Renormalisation group In statistical physics, the renormalisation group method systematically 
implement some form of coarse- graining operation (e.g. rescaling length scales with a factor b: 
L — > — > 4j to expose the character of the large-scale phenomena, in systems where many length 
scales are important. 

Scale invariance A physical system is said to exhibit scale-invariance if it remains unchanged (in 
a statistical sense) by a coarse- graining operation. 

Universality The phenomena whereby many microscopically quite different physical systems ex- 
hibits critical point behaviour, with quantitatively identical features, such as critical exponents. In 
percolation theory e.g. the critical exponents do not depend on the microscopic details of the lattice 
but only on the Euclidean dimension d. 
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